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Abstract. During accretion, a neutron star (NS) is spun up as angular momentum is transported through its liquid surface 
layers. We study the resulting differentially rotating profile, focusing on the impact this has for type I X-ray bursts. The viscous 
heating is found to be negligible, but turbulent mixing can be activated. Mixing has the greatest impact when the buoyancy at 
the compositional discontinuity between accreted matter and ashes is overcome. This occurs preferentially at high accretion 
rates or low spin frequencies and may depend on the ash composition from the previous burst. We then find two new regimes 
of burning. The first is ignition in a layer containing a mixture of heavier elements with recurrence times as short as ~ 5 — 30 
minutes, similar to short recurrence time bursts. When mixing is sufficiently strong, a second regime is found where accreted 
helium mixes deep enough to bum stably, quenching X-ray bursts altogether. The carbon-rich material produced by stable 
helium burning would be important for triggering and fueling superbursts. 
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INTRODUCTION 

As neutron stars (NSs) in low mass X-ray binaries ac- 
crete material from their companions they are expected 
to be spun up by this addition of angular momentum, 
possibly becoming millisecond pulsars This suspi- 
cion has received support by the discovery of accretion 
driven millisecond pulsars jlSll . as well as the millisec- 
ond oscillations seen during type I X-ray bursts 
unstable ignition of the accumulating fuel ifTol 
Angular momentum must be transported into the NS in- 
terior if it is to be spun up, which implies a non-zero, 
albeit small, shear throughout the outer liquid parts of 
the NS. 

Such shearing may lead to viscous heating and chemi- 
cal mixing at depths far below the low density boundary 
layer where the majority of the shearing occurs. This has 
motivated us to assess the importance of angular momen- 
tum transport in NS surface layers. In the following we 
summarize the results of our study It 1 1|] . 



ANGULAR MOMENTUM TRANSPORT 

Material accreted at a rate M reaches the NS sur- 
face with a nearly Keplerian spin frequency of i^K = 
{GM/R^f^ = 1.4 X 10"* s-^ m\'];R^^'^, where M1.4 = 
M/\AMq and R(,=R/\Qf'' cm, which has a kinetic en- 
ergy per nucleon of k, 200 MeV nuc^^ The majority of 



this energy is dissipated in a boundary layer of thickness 
^^BL ^ R and never reaches far into the NS surface 
Nevertheless, angular momentum is added at a rate of 
MR^Q.ji, so that a torque of this magnitude must be com- 
municated into the NS. This implies a non-zero shear rate 
in the interior liquid layers, down to the solid crust. In the 
present work we are interested in the shear at the depths 
where X-ray bursts ignite, near p w 10^ g cm^^, which 
is well below the boundary layer. 

The viscous timescale, fyisc = H^/v, where H is the 
pressure scaleheight and v is the viscosity, is always 
much shorter than the accretion timescale, face = y/in, 
where y is the column depth and rii is the mass accretion 
rate per unit area. Therefore we assume that angular 
momentum is transported in steady-state. In the limit 
i2 < i^K and H -^R this results in Q 



4nR^pvqQ. = MR^Q.K, 



(1) 



where q = dlogD./dlogr is the shear. This equation 
shows that ^ > when angular momentum is transported 
inward. In general, we find that the shearing is rather 
small ^ 1) at the depths of interest. Nevertheless q is 
large enough to activate instabilities that help to transport 
angular momentum, as well as mix material. 



TURBULENT MIXING 

We use equation ([T]i to evaluate the shear for various 
turbulent viscosities and find that purely hydrodynamic 



instabilities ^ are insufficient to prevent strong shear- 
ing of magnetic fields. This leads to generation of the 
Tayler-Spruit dynamo where toroidal field growth 
is balanced by Tayler instabilities to create a steady-state 
magnetic field, which provides a torque on the shearing 
layers. The result is nearly uniform rotation and little vis- 
cous heating, too little to affect either X-ray bursts or 
superbursts, thermonuclear ignition of ashes from previ- 
ous X-ray bursts Ii5j JJJ. Turbulent mixing is found to be 
non-negligible and in some cases may mix fresh material 
with the ashes of previous bursts, therefore we focus on 
this for the majority of our study. 



The Buoyancy Barrier 

A key point for whether mixing can occur is whether 
the strong buoyancy due to the larger density of the ashes 
below can be overcome. Analytic analysis shows that for 
this to occur the accretion rate must exceed 

«W,«5xlO ,„Edd«rs(^^j (^^j, (2) 

where mndd = 1-5 x 10^ gcm^^s^^7?g\ is the local 
Eddington limit for helium-rich accretion, ajs is a di- 
mensionless constant (of order unity) that adjusts the 
strength of mixing, and Aln/x is the fractional change 
in mean molecular weight across the boundary (« 0.44 
for helium-rich material above an iron-rich ocean). This 
critical m depends on the composition of the ashes from 
previous bursts (which sets A In /i) so that the strength of 
mixing could very well change from burst to burst. 



Mixed Unstable Ignition 

We next assess the affect on the burst properties when 
mixing occurs. Figure [T] shows an example calculation 
of material accumulating and mixing on the NS surface 
causing premature ignition of the helium, for the param- 
eters m = O.lmedd and Q. = O.li^K. The shear profile is 
evaluated using equation ([T), and we assume that mix- 
ing occurs completely down to a depth y^ix at which 
face — tmix (the right end of each thin solid line), where 
fmix — H^/D and D is the turbulent mixing diffusivity 
of the Tayler-Spruit dynamo [12]. This is significantly 
deeper than the amount of material accreted, y^cc (shown 
by the filled circles). This means that helium is signif- 
icantly depleted at the time of ignition by an amount 
Ymix — yacc/vmix, where we use a helium-rich composi- 
tion for the accreted material. 

In Figure |2] we plot the recurrence times for mixed 
bursts (once again assuming that eq. [2] is satisfied). 
These timescales are similar to those found for short 
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FIGURE 1. The four panels show how the fully mixed ac- 
cumulating layer evolves in time until it reaches conditions 
necessary for unstable ignition. In each panel, the column of 
helium that has been accreted, y^co is denoted by a filled circle. 
Mixing takes place down to the column reached by the thin 
solid line, ymix. The mixed helium fraction, K^ix = .Vacc/ymix is 
displayed in the upper left-hand comer of each panel. Helium 
ignition curves are shown as a thick dashed line. 

recurrence time bursts To explain the energetics of 
these bursts (i.e. their low a-value |2]) still requires in- 
complete burning from the previous burst. This is in fact 
preferential for strong mixing since incomplete burning 
also leads to smaller compositional gradients. If the mix- 
ing is too strong we instead find stable helium burning 
(noted in Fig.|2]l. This is analytically estimated to occur 
at 

/ £2 \ ^-^^ 

'«crit,2 ~ '«Edd Ors"'^^ V £2^ ) ' 
Above this rh we expect X-ray bursts to cease altogether. 

CONCLUSIONS 

Turbulent mixing is sufficiently large to have important 
consequences for X-ray bursts. We constructed simple 
models, both analytic and numerical, to explore mix- 
ing for pure helium accretion [11]. From these mod- 
els we can make a few conclusions that are likely gen- 
eral enough to apply to most viscous mechanisms. As a 
guide, we show the different burning regimes we find in 
Figure[3] These can be summarized as follows: 

• Mixing has trouble overcoming the buoyancy bar- 
rier at chemical discontinuities. Incomplete burning 
results in small compositional gradients so that mix- 
ing is important in subsequent bursts. 
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FIGURE 2. The recurrence time for mixed-ignition models 
as a function of ih. The symbols denoted different spins, as 
shown in the key. Models that are at sufficiently high m or low 
fl do not ignite unstably, and thus are not plotted. 
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FIGURE 3. The three regimes of burning found for mod- 
els including mixing. The boundaries between each regime 
are shaded to emphasize uncertainty in the strength of mixing 
(cuts = 0.7 — 1.5). The light shaded region divides between 
stable and unstable mixed burning (mcrit,2> 6q (3)). The heavy 
shaded region corresponds to the compositional barrier, which 
prevents mixing («Jcrit.l> (U)- Its position can move sig- 
nificantly depending on the composition of material from the 
previous burst. 



Mixing is strongest at large m (when angular mo- 
mentum is being added at greater rates) and small 
Q. (which gives a larger relative angular momentum 
between the NS and accreted material). 
Mixing can lead to two new effects. First, the layer 
may ignite, but in a mixed environment with a 



shorter recurrence time. Second, for strong mix- 
ing accreted helium can mix and burn in steady- 
state, quenching bursts. Both regimes have observed 
analogs, namely the short recuiTence time bursts |2!] 
and the stabilization of bursting seen near one-tenth 
the Eddington rate |01 ■ 

As a future extension of this work, realistic ash compo- 
sitions from previous bursts should be incorporated into 
simulations, so that one can test whether mixing is a suit- 
able explanation for short recurrence time bursts. Mixed 
ignition can be studied using current numerical experi- 
ments 1 16] by just artificially accreting fuel with a mixed 
composition of y,mx w 0.1 — 0.6. These calculations are 
simplified by our conclusion that shearing and heating 
can be ignored, at least for initial studies. Another im- 
portant extension would be to use different compositions 
for the accreted material. The current models assume it is 
pure-helium. Including hydrogen would allow compari- 
son to a wider range of systems. Finally, the carbon-rich 
ashes from periodic stable burning may be important for 
explaining the recurrence times of superbursts ll5llT4ll6l]. 
which also deserves further investigation. 
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